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1.0  INTRODUCTION 


In  most  practical  applications  of  ultrasonics  we  are  faced  with 
the  problem  of  understanding  the  response  of  an  interface  (usually  a 
water-solid  interface)  to  an  Incident  longitudinal  wave  in  water. 

If  deep  targets  lying  many  wavelengths  within  the  solid  are  to  be 
investigated,  the  amount  of  energy  penetrating  the  Interface  is  often 
of  Interest,  whereas  if  the  region  to  be  investigated  is  near  the 
interface  (within  roughly  one  shear  wavelength) ,  the  detailed  response 
of  the  Interface  must  be  understood. 

The  work  reported  here  has  focused  on  the  latter  problem  and,  in 
particular,  on  the  practical  application  of  the  so-called  Rayleigh 
critical  angle  phenomenon  [1-13].  By  properly  adjusting  the  angle  of 
incidence  of  an  incident  longitudinal  wave  in  water,  it  is  possible 
to  generate  a  Rayleigh-type  interface  wave  that  propagates  along  a 
water-solid  interface.  Such  an  Interface  wave  penetrates  one  shear 
wavelength  into  the  solid  [11]  and  as  it  propagates  along  the  interface 
it  radiates  energy  back  into  the  water  [14].  For  this  reason,  such 
waves  are  often  called  "leaky  surface  waves." 

The  radiated  component  has  considerable  practical  significance 
since  it  is  modulated  by  subsurface  information  (information  arising 
in  the  solid  within  the  penetration  depth  of  the  interface  wave  of 
roughly  one  shear  wavelength).  Proper  examination  and  interpretation 
of  this  radiated  component,  which  is  always  unavoidably  mixed  with  a 
specularly  reflected  component,  can  be  expected  to  provide  useful 
information  on  the  character  of  the  solid  material  near  the  interface. 

It  is  for  this  reason  that  the  Rayleigh  critical  angle  phenomenon 
holds  promise  as  a  tool  for  materials  evaluation  and  nondestructive 
testing.  Near  surface  defects,  flaws,  cracks.  Inclusions,  etc.  are 
easily  detected,  characterized,  and  also  Imaged  by  observing  the 
influence  of  such  anomalies  on  the  radiated  wave  component. 

Although  such  "interface  problems"  have  been  the  subject  of 
numerous  experimental  and  theoretical  studies,  the  Inherent  complexity 


of  the  phenomenon,  when  all  relevant  physical  features  are  Included, 
makes  complete  interpretation  of  experiments  very  difficult. 

In  this  report  we  will  describe  a  new  type  of  experimental 
apparatus  [12,  13]  which  eliminates  at  least  one  of  the  data  inter¬ 
pretation  problems  by  making  possible  local  rather  than  global 
measurements  of  the  Interface  response. 

The  present  report  begins  by  reviewing  appropriate  theory  and 
describing  some  of  the  problems  involved  in  interpreting  Interface 
data  when  the  Interface  is  treated  in  a  realistic  physical  fashion. 

The  nontrivial  Inverse  problem,  namely  going  from  experimental  data 
to  the  elastic  properties  of  the  solid,  is  then  discussed.  The  effects 
of  the  focused  sound  source  used  to  excite  Interface  waves  locally, 
including  nonlinear  effects,  are  also  described. 

We  then  turn  to  a  description  of  the  new  apparatus  and  to 
experiments  on  materials  ranging  from  those  with  isotropic  elastic 
properties  to  those  with  anisotropic  elastic  properties. 

'Images'  of  the  phase  and  amplitude  of  the  measured  local 
reflection  coefficient  (corresponding  to  near  surface  defects)  made 
by  scanning  the  focal  region  of  the  acoustic  lens  over  the  interface 
are  then  described.  The  report  closes  with  a  discussion  of  the 
prospects  and  problems  associated  with  further  practical  development. 


2.0  THEORETICAL  BACKGROUND 


The  correct  (complete)  classical  description  of  the  mechanical 
response  of  a  water-solid  interface  would  be  a  very  complex  affair. 
Such  a  description  would  correctly  treat  all  nonlinear  effects  (in 
the  water  and  the  solid)  which  includes  the  detailed  description  of 
the  transfer  of  linear  and  angular  momentum  at  all  angles  of  incidence 
and  at  all  amplitudes  however  small.  It  would  correctly  treat 
anisotropy,  nonuniformity,  inhomogenity,  and  frictional  attenuation, 
among  other  effects. 

Owing  to  such  complexities,  the  materials  on  both  sides  of  the 
interface  are  often  treated  as  if  they  were  linear,  isotropic, 
uniform  and  homogenous,  and  in  many  practical  applications,  these 
assumptions  are  valid  approximations  to  reality. 

While  the  simplification  to  a  linear  theory,  in  particular, 
leads  to  an  Incomplete  and  sometimes  contradictory  formalism,  the 
practical  significance  of  this  situation  for  many  important  appli¬ 
cations  is  not  great.  Except  for  the  nonlinear  effects  seen  in  the 
water  path,  our  experiments  are  well  described  by  linear  equations 
of  motion  provided  anisotropy  is  also  included  where  appropriate. 


2.1  Anisotropy 

All  crystalline  solids  and  many  alloys  are  anisotropic  [16, 
S3-26].  This  leads  to  equations  of  motion  involving  the  tensor 
properties  of  a  material.  If  these  anisotropic  materials  are  assumed 
to  be  described  approximately  by  linear  equations  of  motion  (all 
anharmonlc  effects  dropped),  the  complete  classical  strain  tensor, 
which  Includes  nonllnearltles ,  namely 


3X,  3X 


and  the  effects  of  anisotropy  are  evidently  contained  in  the  elastic! tv 

tensor  X  „  . 

ikHni 

This  equation  represents  three  homogeneous  equations  in  the  un¬ 
known  particle  displacement  U  =  U  ,  U.  =  U  ,  and  U  =  U  .  Nontrivial 

*  X  2  y  3  Z 

solution  requires  that  the  determinant  of  the  coefficients  vanish 


l-'mrn  '‘k  •'i  -  =  0  . 


(8) 


This  last  equation  is  seen  to  be  a  cubic  in  w’  since  i,  k,  P, 

=  1»  2,  3.  Each  of  the  three  roots  of  this  equation  provide  a 
functional  relation  between  frequency  w  and  wavenumber  k.  Hence, 
there  are,  in  general,  three  kinds  of  independent  motion  and  three 
characteristic  velocities  in  an  unbounded  anisotropic  linear  solid. 

Two  roots  of  Eq.  (8)  correspond  to  bulk  quasi-transverse  waves  and 
one  to  a  bulk  quasi-longitudinal  wave. 

The  existence  of  multiple  bulk  waves  in  an  anisotropic  solid 
greatly  complicates  the  interface  problem.  In  particular,  it  is 
found  that  two  Rayleigh  type  interface  waves  are  generally  produced, 
whereas  only  one  Rayleigh  type  wave  is  produced  at  a  water  Isotropic 
solid  interface  [27-32].  Moreover,  energy  flux  directions  generally 
deviate  from  wavevector  directions  in  anisotropic  solids,  making  simple 
data  analysis  based  on  Snell's  law  (as  applied  to  water  isotropic  solid 
interface)  generally  invalid. 

2 . 2  Interface  Waves  in  Linear  Isotropic  Materials 

The  foregoing  description  of  a  linear  anisotropic  material  is  in 
basic  agreement  with  experience  for  many  applications.  This  general 
formulation  includes  the  special  case  of  a  linear  isotropic  solid 
where  the  equation  of  motion  in  vector  notation  is  approximated  by  [16] 

DU  =  pV^U  +  (X  +  p)V(y*U)  (9) 


-5- 


and  where  X  and  u  are  the  two  Lame  coefficients  which  can  also  be 
written  in  terms  of  the  Young's  modulus  E  and  Poisson's  ratio  o  as 


""  2(1  +  o)(l  -  2o) 

=  E 
2(1  +  o) 

In  terms  of  the  second  order  elastic  constants  ,  we  also  have 

Cii  =  C22  =  C33  =  (X  +  2u) 

C12  =  C21  =  Cl  3  =  C31  =  C23  =  C32  =  X 

Ci,4  =  Css  =  Ces  =  U  • 

In  such  a  linear  isotropic  material,  it  is  possible  to  think 

of  the  particle  displacement  U  as  being  decomposed  into  just  two 
kinds  of  independent  motions  because  of  linearity,  namely 

U  =  V(})  -f  vx? 


where  4)  and  'V  are,  respectively,  the  scalar  and  vector  "potentials' 
describing  the  displacement  U.  Thus,  U  consists  of  a  dilational 
part  Vet)  (characterized  by  a  volume  change)  and  a  rotational  part 
^X'V  (volume  invariant).  Substitution  of  U  into  Eq.  (9)  gives  two 
equations  of  motion,  one  for  4:  and  one  for  'F,  namely 


=  0 


( 


where  C^^  is  the  bulk  longitudinal  wave  velocity 
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where  is  the  bulk  transverse  wave  velocity 
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2p(l  +  0)  • 
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(19) 


2. 3  "Rayleigh"  Critical  Angle  Measurements 

In  addition  to  two  bulk  wave  solutions  in  an  isotropic  solid, 

there  are  surface  wave  solutions  on  a  free  surface  (no  water  loading) 

which  propagate  at  the  velocity  C_  given  by 

K 

Cj^  =  c^n  (20) 

where  n  is  the  single  root  of  the  Rayleigh  characteristic  equation  [11] 
r,®  -  8n“  +  8n'  ^3  -  2  -  16^1  -  0  (21) 

having  the  value  Ini  <  1.  This  is  given  approximately  by  [11] 

=  .87  +  1.12a 

1  +  0  .  (22) 

Here,  C,j,  and  are,  respectively,  the  transverse  and  longitudinal 
velocities  and  a  is  Poisson's  ratio.  These  waves  are  the  so-called 
Rayleigh  waves  that  propagate  on  a  vacuum-solid  interface.  They  are 
characterized  by  a  rapid  reduction  in  amplitude  as  a  function  of 
distance  from  the  boundary  and  penetrate  roughly  one  wavelength  below 
the  solid  surface. 

Even  with  liquid  loading,  the  surface  waves  propagate  along  the 
liquid-solid  interface  with  essentially  the  same  velocity  C„.  However, 
when  liquid  loading  is  present,  a  new  and  very  useful  feature  enters 
as  illustrated  in  Fig.  1.  The  propagating  surface  waves  "leak"  back 
into  the  water,  combine,  and  interfere  with  the  specularly  reflected 
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Figure  1.  Focused  acoustic  source  and  point  receiver  combination  used  in 
Mperiments.  The  source  has  an  aperture  of  2.54  cm  and  a  10  cm  focal  length 
^e  transducer  had  a  5  MHz  center  frequency  and  is  typically  driven  at  10^ 
volts  peak-to-peak  (at  various  frequencies  1-10  MHz).  Amplitudes  recorded 

typically  in  the  1  volt  range  away  from  the 
critical  angle  while  they  are  typically  in  the  millivolt  range  at  the 
critical  angle. 


Various  Investigators  have  studied  the  resultant  nonspecular 
reflection  field  obtaining  expressions  for  the  reflection  coefficient 
under  various  conditions  [28],  [5].  In  general,  the  reflection 
coefficient  for  a  spatially  unbounded  beam  may  be  written  following 
Richardson  [5]  as  (linear  Isotropic  solid) 


cos^2y22  +  (c??/ci 7)^sin2Yi 2Sin2Y22  +  pi ci i sinyi ? /r?ci zsinyi i 


cos22y22  +  (c22/ci2)^sln2Yi 2sln2Y22  -  PiCiisinYi2/p2Ci2sinYii 


where  cj? 


C22 
Cl  1 
■  Yii 

Yl  2 
Y22 
Pi 
P2 


velocity  of  longitudinal  waves  in  the  solid, 

velocity  of  shear  waves  in  the  solid, 

velocity  of  longitudinal  waves  in  the  liquid, 

angle  of  incidence  of  the  longitudinal  wave  in  the 
liquid , 

angle  of  refracted  longitudinal  wave  in  the  solid, 
angle  of  refracted  shear  wave  in  the  solid 
density  of  liquid, 
density  of  solid. 


Angles  are  measured  in  the  counterclockwise  direction  from  the  plane 
of  the  Interface.  The  velocities  are  the  complex  quantities 


c  =  v(l  +  iva/cj)  , 


(24) 


V  =  normal  velocity, 
a  =  attenuation, 

CO  =  radian  frequency. 

The  refracted  angles  are  found  by  the  generalized  Snell's  law. 


cosY^j 

‘'ij 


h  . 


(25) 


The  reflection  coefficient  exhibits  a  sharp  dip  and  a  phase 

shift  at  the  so-called  Rayleigh  critical  angle  6  where  (C, ,  =  V  = 

K  W 

longitudinal  wave  velocity  in  water) 

V 

Sine„  =  .  (26) 


The  above  discussion  applies  to  isotropic  materials  such  as,  crystals 
with  small  anisotropy  parameters  or  glasses  and  polycrystalline  alloys 
with  no  appreciable  residual  stresses  or  grain  alignments.  When 
anisotropies  are  present,  the  problem  is  much  more  complex  since 
Eq.  (26)  breaks  down  and  the  phenomenon  of  multiple  surface  wave 
critical  angles  occurs. 


3.0  QUALITATIVE  ANALYSIS  OF  THE  EFFECT  OF  CONVERGING  BEAMS 


Figures  1  and  2  illustrate  the  basic  experimental  arrangement 
and  the  associated  electronic  apparatus  including  an  IBM  personal 
computer  which  controls  the  system.  The  heart  of  the  system  is  the 
goniometer  which  is  provided  with  a  computer  controlled  stepper  motor 
for  changing  the  angle  of  incidence  (6)  of  the  acoustic  lens  axis. 

Data  analysis  and  display  functions  are  also  carried  out  by  the  computer 
Data  is  stored  on  disk,  displayed  on  a  CRT  and/or  printed  out  on  paper. 
The  system  possesses  a  number  of  Interesting  capabilities  including  the 
ability  to  produce  scanned  Images  of  the  Interface  properties.  We  may 
also  perform  spectral  analyses  of  the  reflected  and  radiated  signals. 

The  convergent  beam  of  sound  is  different  than  other  measurement 
methods,  however,  and  the  question  can  certainly  be  raised  as  to  whether 
our  measurements  are  comparable  to  those  made  by  different  means.  In 
this  section  we  provide  an  initial  analysis  of  our  experiment  based  upon 
well  known  optics  and  acoustics  theory. 


Fig.  2.  Experimental  Apparatus;  IBM  computer  and  display  (center), 
goniometer  and  water  bath  (right),  and  electronics  (left). 
Temperature  for  all  measurements  was  maintained  at  20°C  by 
means  of  the  water  bath  and  the  insulated  tank  at  the 
right  in  this  figure. 


3.1  Plane  of  Incidence 


I 


i 

■ 
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In  order  to  obtain  high  resolution  critical  angle  images  as 
described  in  our  initial  proposal,  we  chose  the  geometry  shown  in 
Fig.  3.  The  transmitter  crystal  is  followed  by  a  lens  which  focuses 
the  sound  to  a  point  at  the  center  of  rotation  of  the  goniometer.  The 
sample  is  placed  with  its  surface  also  at  the  center  of  rotation,  and 
its  surface  normal  in  the  plane  of  the  goniometer.  The  receiving 
transducer  consists  of  a  very  small  element  placed  on  the  other  arm 
of  the  goniometer.  The  central  ray  of  the  transmitter,  and  the 
receiver  fulfill  the  specular  reflecting  geometry. 

The  experiment  is  performed  by  energizing  the  transmitter  with  a 
continuous  sine  wave  at  the  approorlate  frequency  and  monitoring  the 
receiver  with  a  vector  volt  meter.  The  sample  is  placed  in  the  x-y 
plane  and  the  goniometer  carrying  the  transducers  is  moved  from  a  small 
angle  to  a  large  one  or  vice  versa.  In  this  process,  the  transducers 
pass  through  all  three  critical  angles,  0^^,  6^,  6j^  (longitudinal, 
transverse,  and  Rayleigh,  respectively).  The  complex  reflectivity  as 
a  function  of  angle  is  the  recorded  variable. 

In  the  past,  this  type  of  experiment  has  been  performed  with 

large  transducers  and  plane  waves.  Therefore,  each  of  the  critical 

angles  mentioned  above  was  interrogated  separately,  except  perhaps  for 

6  and  6  which  lie  quite  close  together  for  some  materials.  In  our 
1  R 

geometry,  the  focused  cone  can  sometimes  contain  all  three  critical 
angles.  Therefore,  in  addition  to  the  Rayleigh  surface  wave,  there  may 
be  lateral  waves  which  then  radiate  to  form  head  waves  that  will  be 
seen  by  the  receiver.  For  our  usual  geometry,  a  25  ram  transmitter  is 
focused  at  100  mm  resulting  in  a  converging  cone  of  14.25*  Included 
angle.  In  aluminum,  for  example,  0^^  *  13®  26',  0^  =  29®  20'  and 
0j^  =  31®  36'.  Hence,  0^^  and  9^  could  not  be  encompassed  simultaneously 
although  a  slightly  larger  transducer  could  do  this.  0^  and  6j^ 
are,  however,  often  energized  together.  Thus,  a  Rayleigh  and 


lateral  wave  will  often  be  present  simultaneously.  Both  of  these 
waves  will  generate  head  waves  in  the  water  which  will  be  seen  by 
the  receiver. 

3.2  Plane  of  Solid 


Returning  to  Fig.  3,  we  see  that  surface  waves  are  generated 
over  a  wedge  of  angles  emanating  from  the  focal  point  on  the  surface. 

This  wedge  of  surface  waves  reradiates  a  head  wave  back  into  the  water, 
and  therefore,  must  make  a  contribution  to  the  detected  signal. 

The  wedge  of  surface  waves  predicted  by  geometric  theory,  as 
shown  in  Fig.  A,  will  subtend  an  angle  H*'  which  depends  on  the  angle  of 
incidence  of  the  lens  axis  0,  and  the  angle  subtended  by  the  lens 
aperture  from  the  lens  focus  'f. 

Thus,  the  wedge  of  surface  waves  always  has  an  opening  angle  'f' 
which  is  somewhat  larger  than  the  opening  angle  of  the  lens  T. 

From  the  foregoing,  it  is  clear  that  we  cannot  avoid  contributions 
at  the  detector  from  waves  radiated  by  the  wedge  of  surface  waves. 
However,  it  is  not  evident  that  this  will  lead  to  results  differing 
significantly  from  those  expected  for  incident  plane  waves. 

From  a  geometrical  standpoint,  we  note  that  a  wedge  of  surface 
waves  will  produce  a  conically  shaped  diverging  wave  front,  as  shown 
in  Fig.  A,  of  radiated  waves  which  impinge  on  the  point  detector. 

These  conical  surfaces  will  be  characterized  by  a  cone  with  opening  angle 
'f  and  an  axis  parallel  to  the  axis  of  the  lens  in  the  plane  of  incidence. 

The  specularly  reflected  incident  rays  will  produce  a  diverging 
spherical  wave  front  which  also  impinges  on  the  same  detector.  The 
phase  difference  between  the  conical  (nonspecular  reflected)  and  spher¬ 
ical  (specularly  reflected)  wave  fronts  in  the  plane  of  incidence  should 
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LOCUS  OF  SHEAR  CRITICAL  ANGLE 


be  the  same  as  the  phase  difference  between  an  incident  specularly 
reflected  plane  wave  and  its  nonspecularly  reflected  counterpart 
which  is  also  plane.  Away  from  the  plane  of  incidence,  of  course, 
this  statement  no  longer  holds.  Thus,  if  a  very  small  point  detector 
is  used,  the  effect  of  off-axis  rays  should  be  greatly  reduced  if  not 
eliminated  entirely. 

For  example,  in  experiments  on  quartz  and  copper,  we  obtain  results 
that  are  very  similar  to  those  of  other  investigators  who  assumed 
incident  plane  wave  fronts.  Thus,  we  conclude  that  even  for  aniso¬ 
tropic  materials  (which  will  distort  the  radiated  conical  wave  fronts) 
off-axis  excitation  will  probably  not  have  a  major  effect  on  measure¬ 
ments  involving  a  point  detector. 
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4.0  EXPERIMENTAL  RESULTS 


The  results  of  measurements  performed  on  a  4340  steel  specimen 
are  illustrated  in  Fig.  5.  See  Fig.  6  for  a  description  of  the  angles 
involved.  From  Fig.  5  it  might  at  first  appear  that  the  entire  specimen 
is  isotropic.  That  is,  since  the  critical  angle  is  close  to  0  -  31.5° 
for  all  orientations  (J),  there  would  appear  to  be  only  two  elastic 
constants  (X  and  u  or  E  and  o),  and  the  equation  of  motion  would  be 

pU  =  yV^U  +  (X  +  y)y(V*U)  (28) 

admitting  only  Rayleigh  type  surface  wave  solution  [11].  Thus, 
only  one  Rayleigh  type  critical  angle  would  be  expected.  Note, 
however,  that  these  measurements  characterize  a  region  on  the  liquid- 
solid  interface  which  measures  roughly  only  one  wavelength  in  water 
(X  =  .04  cm)  in  diameter. 

We  have  found  other  regions  on  this  same  specimen  (not  illustrated) 

that  are  anisotropic  and  exhibit  more  than  one  critical  angle.  This 

observation  demonstrates  how  the  focused  lens  source  and  point  receiver 

combination  allow  us  to  "map"  or  "image"  the  local  elastic  properties 

at  different  points  on  a  surface,  in  principle.  A  true  isotropic 

material  namely,  optical  glass,  also  shows  a  single  critical  angle  at 

9  =  28.5°  as  illustrated  in  Fig.  7.  Unlike  the  4340  steel  specimen, 

cr 

however,  there  are  no  regions  of  appreciable  anisotropy  in  the  glass 
specimen  so  that  one  critical  angle  characterizes  the  entire  specimen. 

Figures  8  and  9  Illustrate  data  taken  on  a  single  crystal 
of  copper.  This  Illustrates  the  ability  of  the  measurement  system 
to  orient  such  a  crystal.  Since  copper  is  cubic  (see  Fig.  10), 
rotation  about  the  (100)  axis  would  be  expected  to  produce  critical 
angles  showing  the  90°  periodicity  (with  4>)  observed  in  Fig.  8. 

Thus,  the  axis  of  rotation  was  in  fact  the  (100)  axis  and  the 
copper  crystal  has  been  oriented  using  critical  angle  measurements. 


Critical  angles  (6)  vs  (J)  for  a  43A0  steel  sample.  The  frequency 
was  4.1  MHz,  temperature  20®C,  and  source  oscillator  amplitude 
6.5  volts.  The  material  appears  to  be  approximately  isotropic 
in  the  small  region  of  measurement  (the  lens  focus). 


Rotation  of  Samples.  Samples  are  arbitrarily  marked  with 
a  spot  as  shown.  The  spot  is  located  along  the  positive 
x-axis  which  initially  lies  in  the  incident  plane  (S-O-R) 
of  the  source  and  receiver.  The  normal  to  the  specimen 
surface  lies  along  the  z-axis  for  all  azimuthal  angles  4>. 
We  arbitrarily  label  the  x  direction  (J)  =  0®  and  turn  the 
specimen  counterclockwise. 


CRITICAL  ANGLE  0 


cr 


Critical  angle  for  glass.  The  figure  shows 
the  critical  angle  as  a  function  of  the 
azimuthal  angle  ip  (angle  of  rotation  of  the 
specimen) . 


90* 


80* 


Raylelgh-type  (surface  wave)  critical  angles  for  a 
copper  single  crystal  rotated  about  the  (100)  axis. 
The  frequency  was  3.8  MHz  and  the  two  critical  angles 
occur  at  0c  -  41*  and  0c  =  47*.  The  azimuthal  angle 
^  (angle  of  rotation  about  the  100  axis)  varied  from 
0  to  360*  and  data  was  taken  in  steps  of  A(t)  «  S*. 

Pure  mode  directions  (001,  310,  and  010)  are  marked 
on  the  figure. 
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Fig.  9.  Amplitude  R(-3)  and  phase  v(h)  vs  P  for  different  azimuthal  angl 
copper  single  crystal  (angle  of  rotation  ‘  about  th*-  100  axis). 
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(b)  (100)  axis 


(c)  (110)  axis 


(a)  (111)  axis 

(111)  plane  edge  on 


g.  10.  Miller  indices  for  cubic  crystals.  The  axes  are  perpen¬ 
dicular  to  the  shaded  planes  of  the  same  designation. 
That  is,  the  100  axis  is  perpendicular  to  the  100  plane. 
Many  directions  are  qulvalent  in  a  cubic  crystal.  For 
example,  the  100  direction  is  equivalent  to  the  001 
direction  or  the  010  direction  in  so  far  as  bulk  elastic 
properties  are  concerned. 


Not  all  crystals  may  be  oriented,  however,  as  illustrated  by 
the  measurements  made  on  an  aluminum  single  crystal.  The  crystal  of 
aluminum,  having  a  very  low  anisotropy  parameter  [24-26],  looks 
very  nearly  isotropic  as  illustrated  by  the  critical  angle  measure¬ 
ments  of  Fig.  11.  Clearly,  such  crystals  cannot  be  oriented  easily 
using  critical  angle  measurements  alone. 

Finally,  in  Fig.  12  we  illustrate  a  series  of  measurements 
made  on  a  titanium  alloy  with  varying  O2  content.  One  population  of 
samples  (two  samples)  with  different  but  "low"  O2  content  (A  and  B) 
seem  to  have  more  or  less  isotropic  characteristics.  But,  above  a 
certain  O2  content  (samples  C,  D,  and  E)  all  of  the  samples  appear 
to  be  anisotropic  owing  to  the  distinct  elliptical  patterns  exhibited 
by  curves  of  6  versus  <{i  as  in  Fig.  12.  Although  the  data  is  sketchy 
(more  points  on  a  given  sample  should  be  examined  and  averaged,  for 
example),  it  does  appear  that  a  change  in  crystal  structure  has  taken 
place  above  a  threshold  O2  content.  In  particular,  below  this 
threshold  the  anisotropy  parameter  is  small  and  above  this  threshold 
it  is  large.  Although  the  detailed  reasons  for  such  behavior  reside 
in  the  domain  of  metallurgy  (wavelength  exceeds  sizes  of  crystallites 
so  that  a  long  range  correlation  or  alignment  over  many  crystallites, 
which  depends  on  O2  content,  appears  to  take  place  at  specimen  C) ,  it 
is  clear  that  acoustic  critical  angle  techniques  might  be  helpful  in 
characterizing  or  signaling  the  onset  of  such  changes  of  the  crystal¬ 
line  symmetry  character  of  a  specimen. 


fig*  11.  Critical  angles  (0)  vs  ((>  for  an  aluminum  single  crystal.  The 
frequency  was  3.8  MHz,  temperature  20“C,  and  source  oscillator 
amplitude  9  volts.  The  axis  of  rotation  was  coincident  with 
the  (100)  axis  of  this  cubic  crystal. 


This  figure  Illustrates  critical  angle  measurements  on  five 
specimens  of  a  titanium  alloy  with  varying  O2  content.  The  O2 
content  increases  from  A  to  E.  Curves  for  A  and  6  seem  to  be 
more  or  less  isotropic  (circular  pattern)  while  C,D,E  fall  in 
a  distinctly  elliptical  pattern  suggesting  that  a  major  crystal 
line  anisotropy  enters  for  specimens  C,D,E.  This  change  in 
crystal  structure  appears  to  occur  at  a  certain  threshold  O2 
content  characteristic  of  specimen  C. 


5.0  IMAGING  EXPERIMENTS 


Perhaps  the  most  Interesting  application  of  the  Rayleigh  critical 
angle  phenomenon  is  in  imaging  near  surface  mechanical  properties.  As 
mentioned  earlier,  we  are  not  yet  armed  with  computer  algorithms  that 
would  allow  a  complete  working  out  of  local  elastic  constants  from 
measurements.  However,  one  can  indirectly  detect  changes  in  these 
elastic  properties  by  noting  changes  in  the  reflected  amplitude  R(G) 
and  the  phase  1(6)  at  the  point  detector  (see  Fig.  1).  The  local 
changes  over  the  interface  constitute  an  "image"  of  whatever  (flaw, 
discontinuity,  etc.)  it  is  that  is  responsible  for  a  change  in  the 
elastic  properties  of  the  solid  near  this  surface.  As  an  example  of 
what  can  be  done,  consider  the  images  produced  in  Figs.  13  through  16. 

In  each  case,  except  for  Fig.  16,  the  specimen  was  aluminum,  which  as 
we  have  seen  is  very  nearly  isotropic  even  in  the  crystalline  state 
(see  Fig.  11). 

Because  the  reflected  longitudinal  wave  is  modulated  by  the 
interface  wave  and  because  the  interface  wave  extends  roughly  one  shear 
wavelength  (X-f  =  Cj/f)  into  the  solid,  these  images  provide  information 
on  the  solid  within  X-p  of  the  liquid-solid  boundary. 

The  images  in  Figs.  13  through  16  were  made  by  scanning  the  sample 
(using  a  computer  controlled  x-y  scanner  not  shown  in  Fig.  2).  Targets 
scanned  ranged  from  flat-bottomed  holes  in  aluminum  to  a  piece  of  fine 
grained  granite.  Note,  as  in  Fig.  13,  that  whenever  the  boundary 
conditions  are  insufficient  to  generate  a  surface  wave  (9  ^  9j^>  '^j)  or 
the  depth  of  the  discontinuity  is  greater  than  the  surface  wave  penetra¬ 
tion  (Xj)  as  in  Fig.  14,  no  discontinuity  can  be  seen  in  the  images. 

It 

Figure  15  Illustrates  the  result  of  imaging  a  stress-dislocation 
field  resulting  from  a  steel  ball  bearing  being  pressed  into  an  aluminum 
sheet  [28]  -  [30].  Figure  16  illustrates  images  of  a  granite  specimen 
showing  what  appears  to  be  individual  crystal  grains. 
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0c  +  1.4‘ 


©c  +  1.4' 


Fig.  13.  Phase  ’i'(0)  and  amplitude  R{0)  Images  of  a  .6  cm  diameter  flat- 
bottomed  drill  hole  .8  mm  from  the  surface  of  an  aluminum  plate 
.9  cm.  thick.  Images  at  the  top  of  the  figure  were  made  at  the 
critical  angle  for  aluminum  (=31°)  while  those  below  were  made  at 
successively  larger  angles.  Note  that  as  0  increases,  the  conditions 
for  surface  wave  generation  are  eventually  ruined  and  the  images  grow 
steadily  worse.  The  frequency  for  these  experiments  was  3.75  MHz. 

A  lens  and  point  receiver  similar  to  that  shown  in  Fig.  1  was 
employed.  The  hole  Images  are  somewhat  eliptical  In  shape  owing  to 
the  Incorrect  aspect  ratio  of  the  monitor. 
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PHASE  iriAGES  i/»(o) 


AMPLITUDE  IMAGES  R  (e) 


1.86  MHz  1.86  MHz 


Fig*lA.  Phase  'i'(6)  and  amplitude  R(0)  images  of  a  series  of  .  6  cm  diameter 
flat-bottomed  drill  holes  at  varying  depths  (.A,  .8  and  1.6  mm) 
from  the  surface  of  an  aluminum  plate  .9  cm  thick.  All  images  were 
made  at  the  critical  angle  for  aluminum  (=31°)  and  varying  the 
frequency  as  shown.  Note  that  as  the  frequency  (f)  increases,  the 
depth  of  penetration  of  the  surface  waves  (1-r  C^/f)  decreases, 

and  for  a  frequency  of  10  MHz  (not  shown)  the  hole  at  .A  mm  depth  is 
no  longer  visible  for  example. 
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PHASE  IMAGE  >^(6) 


AMPLITUDE  IMAGE  R(e) 


STRESSED  REGION 


Fig.  15.  Phase  '1^(6)  and  amplitude  R(6)  images  of  a  stress-dislocation 

field  in  a  3.2  mm  thick  aluminum  plate  in  which  a  4.8  mm  steel 
ball  bearing  was  pressed  to  a  depth  of  .8  mm.  The  frequency 
chosen  was  3.75  NHz  but  the  critical  angle  (not  recorded)  was 
that  of  the  flaw  itself  and  not  the  aluminum  (=31°). 


V" 


Fig.  16.  Phase  v(6)  and  amplitude  R(0)  images  of  a  fine  grained  granite 
specimen  (1-2  mm  diameter  grains)  at  a  chosen  incident  angle  of 
29.5®.  The  frequency  was  3.75  MHz.  Both  phase  and  amplitude 
images  show  what  appear  to  be  individual  crystal  grains  (3.8  cm 
scan  length). 


The  qualitative  interpretation  of  these  images  is  clear  but  the 
quantitative  interpretation  has  yet  to  be  carried  out  since  it  involves 
very  complex  changes  In  the  boundary  conditions  over  short  distances 
near  the  discontinuities.  In  principle,  such  calculations  can  be 
made  for  the  amplitude  images,  however  there  are  questions  regarding 
the  phase  Images  and  the  Information  on  which  they  are  based. 

Figures  17  and  18  Illustrate  an  interesting  situation  that  is 
sometimes  encountered.  Phase  shifts  are  typically  360®  on  passing 
through  the  critical  angle.  That  is,  the  phase  normally  goes  through 
a  change  where  it  continually  increases  (or  decreases)  but  returns  to 
0®  (360®).  Sometimes,  however,  a  curve  of  the  type  shown  in  Fig.  17 
occurs.  The  phase  advances,  retards,  and  then  advances  again  returning 
to  zero.  This  effect  can  be  seen  on  the  same  sample  and  does  not  seem 
to  be  related  to  electronic  noise  or  other  similar  problems.  We  also 
know  that  these  phase  measurements  are  unaffected  (directly)  by  higher 


a; 

a. 

c. 

E 

< 


0/ 

X 

re 

X 

c~ 


Descriptions 

File  name 
Specimen  name 

Frequency  or  other  comments 

Initial  angle  6 
1°  increments  (0) 

Number  of  angular  positions 

Critical  angle  (R(0)  minimum) 
Phase  range  on  'v(0)  graph 
Spacing  on  V  (6)  graph 
Number  of  measurements  averaged 
at  each  angular  position 


Incident  Angle  0  (Degrees) 


Flg>  17.  Reflected  Amplitude  R(9)  and  Phase  i^(6)  for  a  Cu  Single  Crystal 
(100  Plane)  at  20‘’C  and  a  Frequency  of  3.0  MHz.  These  curves 
(|J'(6))  illustrate  the  phase  "reversal"  effect.  Compare  with 
a  "normal"  type  phase  curve  'I' (6)  such  as  that  illustrated  in 

Fig-  18. 


Fig.  18.  Reflected  Amplitude  R(6)  and  Phase  1(^(6)  for  a  Cobalt  Powder 
Composite  Showing  a  "Normal"  Phase  Curve  lii  (6)  at  20°C  at  a 
Frequency  of  3  MHz.  Note  that  the  phase  shift  for  these 
"normal"  phase  curves  is  typically  360°,  whereas  the  curves 
showing  the  "reversal"  (Fig.  17)  is  less.  The  critical  angle 
for  this  specimen  is  located  at  A2°.  The  choice  of  the 
particular  materials  illustrated  in  Figs.  17  and  18  is  arbitrary. 
Many  different  materials  illustrate  these  effects.  The  same 
specimen  often  exhibits  both  types  of  curves. 
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harmonics  since  our  phase  meter  sees  only  the  fundamental.  However, 
we  cannot  rule  out  the  possibility  that  a  signal  at  the  fundamental 
frequency  produced  by  nonlinear  effects  (Intermodulation  of  harmonics, 
for  example)  is  distorting  the  "normal"  phase  curve  near  the  critical 
angle. 

Regardless  of  the  ultimate  explanation  for  such  effects,  it  is 
clear  that  phase  information  and  especially  phase  images  are  not  to  be 
trusted  unless  and  until  these  problems  are  cleared  up.  This  is  not 
the  case  with  the  amplitude  Images  which  would  be  little  altered  by 
contamination  with  a  small  signal  at  the  fundamental  frequency. 


6.0  NONLINEAR  EFFECTS 


We  have  delayed  a  detailed  discussion  of  nonlinear  effects  until 
the  last,  for  two  reasons.  First,  as  emphasized  in  the  beginning  of 
this  report,  many  features  of  the  critical  angle  phenomenon  (such  as 
the  magnitude  of  various  critical  angles)  are  well  accounted  for,  at 
least  in  principle,  by  linear  theory  supplemented  by  anisotropy. 

Second,  the  effects,  that  we  have  observed  which  can  be  attributed  to 
nonlinearities,  are  not  yet  sufficiently  well  understood  to  place  them 
in  a  position  of  prominence. 

It  is  known  that  both  water  and  a  typical  solid  are  nonlinear, 
and  that  water  is  more  nonlinear  (has  a  shorter  discontinuity  length, 
for  example)  than  the  typical  solid  [17],  [18],  [37].  Thus,  we 
expect  to  see  nonlinear  effects  due  to  the  water,  and  possibly  the 
solid  as  well,  at  the  critical  angle  where  the  amplitude  in  the 
solid  increases  dramatically. 

6.1  Spectral  Measurements 

A  dynamical  physical  system,  such  as  the  water  path-solid  combina¬ 
tion  in  our  experiments,  can  be  characterized  by  its  frequency  response. 
That  is,  the  frequency  spectra  of  the  signals  propagating  through  the 
system  may  be  analyzed  to  determine  the  response  of  the  system  to  various 
Inputs.  Any  anomalies,  such  as  nonlinearities,  will  then  be  characterized 
by  changes  in  such  spectra  as  a  function  of  certain  variables  (source 
voltage,  angle  of  incidence,  etc.).  Although  it  is  very  difficult  to 
obtain  absolute  spectral  information  in  a  complex  electronic  and  physical 
system,  owing  to  many  coupled  transfer  coefficients,  it  is  possible  to 
obtain  relative  spectral  Information  as  certain  quantities,  such  as 
incident  angle  0,  are  varied.  We  have  made  such  relative  spectral 
measurements  on  nonspecularly  reflected  longitudinal  waves  (reflected 
and  radiated  by  leaky  surface  waves)  at  and  near  various  critical  angles. 


This  information  provides  insights  into  the  nature  of  the  scattering 
process  as  a  function  of  the  angle  of  incidence  (6)  of  the  longitu¬ 
dinal  waves  and  other  relevant  variables. 

A  typical  spectral  analysis  of  the  reflected  wave  is  performed 
by  first  finding  a  suitable  frequency  at  which  a  sharp  minimum  of 
R(6)  occurs.  A  suppression  filter  is  then  adjusted  to  cut  out  this 
frequency  while  passing  the  higher  harmonics.  It  is  necessary  to 
suppress  the  fundamental  owing  to  the  fact  that  our  digitizer  has  only 
8-bit  accuracy,  and  the  higher  harmonics  are  small  relative  to  the 
unfiltered  fundamental.  Signals  are  digitized  at  a  20  MHz  rate  and 
then  Fourier  analyzed  to  obtain  Che  resultant  frequency  spectrum. 

An  example  of  such  a  spectrum  is  Illustrated  in  Fig.  19. 

Samples  of  mica,  optical  glass  and  aluminum  (single  crystal) 
were  examined  at  the  same  fundamental  frequency  (4.5  MHz)  at  several 
different  source  voltage  settings  using  a  4.5  MHz  suppression  filter. 

6.1.1  Mica  Standard.  The  question  of  whether  nonlinear  effects  due 
to  the  water  (and  possibly  the  solid)  are  important  in  interpreting 
critical  angle  data  is  a  complex  one,  owing  to  the  fact  that  the 
reflection  coefficient  is  capable  of  drastically  altering  the 
harmonic  content  generated  solely  in  water  (not  to  mention  the  solid). 

Consider  the  observations  made  on  a  thin  (.002  in.)  reflecting 
mica  window.  Figure  20(a)  through(c)  illustrates  data  obtained  on 
the  fundamental  (4.5  MHz)  and  the  second  harmonic  (9.0  MHz)  for 
angles  of  incidence  ranging  from  15  to  45°.  Before  the  data  was 
digitized,  it  was  first  passed  through  a  4.5  MHz  (-30  DB)  suppression 
filter  (which  has  very  little  effect  on  the  higher  harmonics) 
amplified  and  then  ilgitized  at  a  20  MHz  rate  (8  bit  accuracy).  The 
data  was  then  Fourier  analyzed  using  an  FFT  routine  (on  the  IBM 
personal  computer)  and  the  amplitudes  were  then  plotted  as  shown  in 
these  figures.  Data  at  three  different  source  voltage  levels  (10,  8, 

6  volts)  is  illustrated  in  these  figures. 


Fig.  20  (h)  Source  voltage  8V 


The  first  thing  to  be  noted  about  the  data  is  that  the  signals 
for  each  voltage  setting  contain  a  second  harmonic  component.  Second, 
the  amplitude  of  both  the  fundamental  and  the  second  harmonic  (at 
each  voltage  setting)  is  essentially  constant  at  its  average  value  for 
any  angle  of  Incidence  6.  Third,  the  ratio  of  the  second  harmonic 
amplitude  to  the  fundamental  amplitude  is  essentially  independent  of 
the  voltage  setting.  This  last  observation  does  not  imply  that  both 
amplitudes  are  independent  of  the  voltage  setting.  In  fact,  the 
average  amplitude  of  the  fundamental  should  be  taken  to  be  roughly 
proportional  to  the  source  voltage  setting. 

These  observed  facts  lead  to  several  important  conclusions 
regarding  the  nature  of  the  source,  the  effects  of  the  mica  reflector, 
and  the  nature  of  the  nonlinear  effects  in  the  water  path.  The  fact 
that  (Aj/Ai)  is  essentially  independent  of  the  voltage  settings  is 
very  significant.  It  means  that  the  effect  of  the  mica  reflector 
completely  cancels  any  second  harmonic  generated  in  the  water  and 
passes  only  the  initial  second  harmonic  contamination  at  the  source. 

To  see  this,  we  first  note  that  an  increase  in  the  voltage  setting 
(call  it  k)  causes  the  fundamental  amplitude  Ai  to  increase  in  the 
same  proportion  as  the  second  harmonic  amplitude  Az ,  thereby 
maintaining  the  ratio  kAz /kAi  =  Az/Ai  =  constant.  Thus,  the  second 
harmonic  we  see  on  the  mica  reflector  is  only  that  due  to  contamination 
at  the  source.  This  is  not  to  say  that  there  are  no  higher  harmonics 
produced  by  the  source.  In  fact,  these  are  present,  but  as  Breazeale 
has  emphasized,  if  each  frequency  component  is  shifted  by  180°  (for 
any  reason)  when  the  waves  return  to  the  receiver  (source  reflector 
distance  =  reflector  -  receiver  distance),  they  will  have  no  higher 
harmonics  present  other  than  initial  contamination.  Now  we  also 
know  that  the  mica  phase  shifts  all  Incident  waves  by  180°  so  that  any 
observed  component  at  the  receiver  is  merely  the  contamination 
present  at  the  source.  This  explains  in  a  simple  way  why  A1/A2  is 
constant  with  source  voltage.  The  source  voltage  (10,  8,  6  volts) 
definitely  produce  higher  harmonics  due  to  nonlinearities,  but  all  of 
these  are  always  eliminated  by  propagation. 
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6.1.2  Measurements  Along  the  Lens  Axis.  That  a  second  harmonic 
component  due  to  nonlinear  effects  is  actually  present  along  the 
water  path  is  demonstrated  by  the  following  experiments.  The 
results  illustrated  in  Figs.  21  through  23  were  obtained  by  examining 
the  signals  along  the  axis  of  the  acoustic  lens  including  the  focal 
region,  which  in  this  case  is  located  at  4.0  inches  from  the  lens 
(geometric  focus). 

Because  the  beam  is  converging,  the  amplitudes  of  both  the 
fundamental  and  any  second  harmonic  contamination  that  originate  at 
the  source  will  generally  Increase  as  the  focus  is  approached,  and 
likewise  would  be  expected  to  diminish  beyond  the  focus.  Note  that 
the  fundamental  component  has  just  such  a  behavior,  l.e.,  it  is 
essentially  symmetric  about  its  maximum  amplitude.  It  should  be 


noted,  however,  that  the  maximum  amplitude  occurs  at  a  position  some¬ 
what  closer  to  the  lens  than  the  geometric  focus. 

The  second  harmonic  (Figs.  22  has  a  considerably  different 
behavior  than  the  fundamental.  Note  that  it  is  not  symmetric  about 
its  maximum  value.  This  can  only  mean  that  a  nonlinear  process  is 
creating  excess  second  harmonic  components  over  and  above  the  input 
contamination.  In  all  of  these  figures,  it  is  clear  that  the 
amplitudes  A2,  Ai  generally  Increase  with  source  voltage  in  a  way 
that  demonstrates  nonlinear  behavior  in  the  water.  For  example,  if 
the  water  were  purely  linear,  a  source  voltage  factor,  k,  would 
result  in  a  constant  ratio  kAz/kAi  =  Aj/Ai  =  constant.  This  is 
clearly  not  the  case.  As  the  source  voltage  is  increased,  this  ratio 
Increases  at  the  detector  for  essentially  all  distances  from  the 
source  indicating  the  nonlinear  production  of  a  second  harmonic. 

6.1.3  Optical  Glass.  The  position  of  the  critical  angle  can  depend 
on  the  frequency  and  this  fact  can  also  affect  the  behavior  of  the 
second  harmonic  amplitude  relative  to  the  amplitude  of  the  fundamental. 
To  see  this,  consider  the  data  taken  on  optical  glass  as  shown  in 
Fig.  24(a)  through  (e) .  Freauency  (fundamental)  and  amplitude  (in  some 


; 
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I 


Distance  from  Lens  (Inches) 


Fig.  22,  Amplitudes  of  the  second  harmonic  (9  MHz)  at  various  distance  from  the 
4.0  In.  focal  length  lens.  Source  voltages  were  10,  8,  and  6  volts. 


Distance  from  Lens  (Inches) 


Ratio  of  second  harmonic  amplitude  (9  MHz)  to  amplitude  of  the  fundamental 
(4.5  MHz)  at  various  distances  from  the  4.0  in.  focal  length  lens.  Source 
voltages  were  10,  8,  and  6  volts. 


15.00 


OPTICAL  GLASS 


Relative  reflected  amplitudes  of  the  fundamental  and  second  harmonic  (9  MHz)  on  optical  gl, 
filter  used  to  suppress  fundamental)  for  different  source  voltages  and  angles  of  incidence 
note  on  amplitude  variations  In  Fig.  36. 


OPTICAL  GLASS 


cases)  were  chosen  to  be  the  same  as  for  the  data  taken  on  the  mica 
window. 

Data  was  taken  for  source  voltage  settings  of  10,  9,  8,  7, 
and  6  volts.  The  data  at  each  voltage  setting  is  characterized  by 
a  second  harmonic  amplitude  which  goes  to  a  minimum  at  an  angle 
ranging  up  to  4°  less  than  the  critical  angle  characterizing  the 
position  of  the  minimum  of  the  fundamental  amplitude.  Moreover,  the 
postion  of  the  second  harmonic  minimum  amplitude  is  constant  as  a 
function  of  the  source  voltage,  whereas  the  position  of  the  minimum 
of  the  fundamental  amplitude  steadily  shifts  to  lower  values  as  the 
source  voltage  is  reduced  (see  Fig.  25). 

6.1.4  Discussion.  It  has  been  stressed  by  Breazeale  and  others 
that  if  each  component  in  an  incident  sawtooth  waveform  is  shifted 
by  180®  and  allowed  to  propagate  without  attenuation,  it  will 
eventually  "heal"  and  be  converted  back  to  the  fundamental  at  the 
detector,  provided  the  detector  is  at  the  same  distance  from  the 
reflector  as  the  source  is  from  the  reflector. 

Now  even  if  we  ignore  frictional  attenuation  (geometric 
attenuation  is  the  same  for  each  frequency  and  can  safely  be  ignored) 
which  differs  for  each  frequency  component,  we  certainly  cannot 
ignore  the  amplitude  modulation  of  each  component  near  critical  angles 
which  vary.  If  the  critical  angle  is  the  same  for  each  harmonic  and 
the  shape  of  the  amplitude  minima  is  identified,  then  the  amplitude 
modulation  that  occurs  when  the  critical  angle  is  approached  should 
have  little  effect  on  the  "healing"  process  described  above. 

In  conclusion,  it  is  not  possible  at  this  time  to  yay  that  we 
have  seen  effects  due  to  nonlinearities  in  the  solid,  but  there  is  no 
doubt  that  nonlinear  effects  in  the  water  coupled  with  the  particular 
properties  of  the  solid  can  significantly  effect  the  resultant 
harmonic  content  in  the  received  signal.  This  is  clearly  an  area 
needing  further  study. 
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7. 0  CONCLUSIONS 


We  have  performed  numerous  critical  angle  experiments  on  a 
variety  of  samples  ranging  from  single  crystals  to  polycrystalline 
alloys  and  glasses.  Two  basic  conclusions  can  be  drawn  from  this 
work:  1)  simple  linear  equations  of  motion  (which  include  aniso¬ 

tropy)  suffice  for  an  approximate  description  of  the  observations 
(for  example,  the  value  of  various  critical  angles  can  be  calculated 
from  a  knowledge  of  second  order  elastic  constants  alone);  2)  certain 
observations  at  and  near  the  critical  angle,  involving  the  production 
of  anomalous  harmonics,  imply  that  the  nonlinear  characteristics  of 
the  water  (and  possibly  the  solid)  cannot  be  ignored  and  need  further 
study.  Moreover,  the  body  of  knowledge  gained  so  far  shows  that  our 
new  measuring  apparatus,  using  an  acoustic  lens,  has  the  following 
desirable  features:  a)  it  allows  local  measurements  of  the  solid 
properties  to  be  observed  and  b)  except  for  problems  associated  with 
phase  measurements  of  the  reflected  waves,  experimental  results  are 
essentially  consistent  with  the  work  of  others  who  used  a  different 
experimental  approach.  We  conclude  that  our  goal  of  using  a  device 
in  practical  flaw  imaging  work  is  reasonable  and  holds  promise  of 
being  a  quantitative  flaw  imaging  technique. 
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